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Abstract 

We discuss three different globally regular non-topological stationary soliton solutions 
in the theory of a complex scalar field in 3+1 dimensions, so-called Q-balls, Q-vortices and 
Q-walls. The charge, energy and profiles of the corresponding solutions are presented for 
each configuration studied. The numerical investigation of these three types of solutions 
shows different behavior of charge and energy with changing frequency uj for each type. 
We investigate properties of new families of coupled non-topological 2-Q-ball solutions 
obtained within the same model by generalization of the ansatz for the scalar field which 
includes an independent phase. New composite solutions for another known model, which 
describes two Q-balls minimally interacting via a coupling term, are discussed briefly. 

1 Introduction 

Q-balls are stationary localized non-topological soliton solutions of a nonlinear field the- 
ory which carry global U{\) charge pj. In order to evade Derrick's theorem [2], Q- 
ball configurations must be time-dependent. In the simplest case there is a single com- 
plex self-interacting scalar field with an explicitly time-dependent phase and sextic non- 
renormalizable potential. The charge then is directly proportional to the frequency of 
rotation. 

As shown in [3j, it is possible to construct spinning axially symmetric generalization 
of the Q-ball solutions. These solutions correspond to stationary localized configurations 
possessing a finite mass and a finite angular momentum which is quantized, J = nQ, where 
Q is the Noether charge of the solution and n G Z corresponds to the winding around the 
symmetry axis. Possessing even or odd parity, their energy density forms one or more 
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tori Recently the properties of these rotating solutions in the presence of gravity were 
considered in [3 E] . 

An interesting class of non-topological solitons having axial symmetry was discussed in 
[3]. These solutions represent stationary Q- vortices. Other configurations which possess 
planar symmetry has been investigated in In the latter case the solutions are refer to 
as Q- walls. 

The Q-ball solutions arise in various models, however one of the most interesting ex- 
amples is related with the supersymmetric extensions of the standard model with fiat 
directions in their scalar potentials [8]. In such a case the U{1) charge is associated with 
the symmetries of baryon and lepton number conservation, so the corresponding solutions 
are leptonic and baryonic balls. It was suggested [9] that the Q-balls arising in a super- 
symmetric model may play a role in baryogenesis through the Affleck-Dine baryogenesis 
mechanism [TU]. 

Since the supersymmetric models contain several scalar fields, they interact via various 
potentials. This pattern can then lead to the complicated picture of evolution of the Q- 
balls. On the other hand, the system of coupled Q-balls may possess new solutions whose 
properties may be rather exotic. For example previous research has suggested that such a 
system may support the existence of 'twisted' Q-balls which would have a certain similarity 
with twisted loops in the Faddeev-Skyrme model [4j. An intriguing observation is that in 
the case of the supersymmetric models the flat Q-ball potential yields the energy/charge 
relation E ~ Q^^'^ [H |T3] , which precisely matches the topological energy bound for the 
solitons of the Faddeev-Skyrme model [H]. Another interesting similarity between these 
models is that the spectrum of solition solutions in both cases includes not only the fun- 
damental localised sohtons, the Q-ball and the hopflons, respectively, but also extended 
objects like vortices [3l [15] and walls [71 ITS] . 

In this paper some properties of Q-ball type solitons are studied. First, we here consider 
properties of the corresponding solutions to the same simple model with non-renormalisable 
|$|^-potential. We discuss both the usual spherically symmetric Q-balls and extended ob- 
jects, Q-vortices and Q-walls. Another purpose of this paper is to investigate properties 
of the 'twisted' Q-balls using the parametrization suggested in PJ. We argue that these 
solutions may exist only if the constituents are actually two copies of the same conflgura- 
tion. Using another model suggested by Brihaye and Hartmann [20| we present numerical 
arguments for the existence of a class of coupled 2-Q-ball conflgurations which may have 
different geometry. We present both spinning and non-spinning solutions and briefly dis- 
cuss their properties. 

The plan of the paper is as follows. In section II, we review the model and give the 
equations and boundary conditions. Here we discuss the Q-ball, Q-vortex and Q-wall 
solutions for the model with sextic potential, while in section III, we discuss our results 
for two different models with two interacting components. We present our conclusions in 
section IV. 
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2 The model and fundamental solutions 



We consider a theory of a complex scalar field $ in (3+l)-dimensional flat spacetime defined 
by the Lagrangian density 

L = a^$9^$*-[/(|$|) (1) 

with a U{1) invariant scalar potential f/(|$|). 

To secure the existence of non-topological solitons, the vacuum of the model should be 
non-degenerate [TT]. In particular, a potential polynomial in |$p should contain powers of 
$ higher than four [1], which makes the model non-renormalizable. We will consider the 
|$|6 potential [3112112] 

f/(|$|) = a|$|6 + 6|<l>|^ + c|$p, (2) 

where a = 1, 6 = -2, c = 1.1, i.e. f/(|$|) = |$p[(l- |$p)^ + 0.1]. This is a typical potential 
in field theories that can contain extended objects of the Q-ball type. So, the mass of the 
scalar excitation is = |f/"(0) = c. 
The corresponding field equation is 

and the global U{1) Noether charge is 

Q = i [ £x{^^* - $*$) (4) 



Depending on the boundary conditions, different types of solutions can be found. The 
well known Q-ball solution corresponds to a minimum of the energy functional at a fixed 
charge p[]. Assuming that $ depends on time harmonically 

$ = 0(r)e^^* , (5) 

where u is the internal rotation frequency and 0(r) is time independent real function of 
coordinates r, we get the total energy functional 

T^ = E = j d'x + {^<Pf + U{<P)) , (6) 

which effectively describes the motion in the potential V^(0) = — U{(j)). The usual 
restriction on the rotation frequency is that the potential U{(f)) should have its absolute 
minimum U{0) = and the Q-ball solution must decay at spatial infinity. This implies [1] 

Cu'i < LO^ < Lol , (7) 

where with our choice of the parameters, = 0.1 and u"^ = m? = 1.1. 

The simplest spherically symmetric Q-ball placed at the origin is a stationary solution 
of the form $ = /(r)e*'^* [Ij, where the real amplitude /(r) satisfies the equation 

^^ + ^^ + c.V = -— . (8) 

dr"^ r dr 2 df 
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The regular solution of this equation asymptotically decays as /(r) ~ re 
charge of the fundamental Q-ball depends on uj as 



The 



(9) 



As — 7- cj^ the Q-ball expands and its charge diverges. In the opposite limit — ?■ u;^ = 
w? the Q-ball also expands. The maximal value of the function /(O) at the origin then 
tends to zero although the integrated charge/energy density diverges (See Fig. [U (a)). 

Another stationary solution of the model ([1]) is the vortex configuration discussed in 
[3]. Assuming the field has cylindrical symmetry, this leads us to the parametrization 



$ = g{p)e^'^\ where p = a/x^ + 



Then the field equation becomes 



£g Idg 2 

dfP pdp ^ 2 dg 



(10) 



with the boundary conditions g'{0) = and g{oo) = 0. So, the field g{p) has a maximum 
at p = and decreases monotonically towards zero at infinity. 






Figure 1: Stationary non-spinning non-topological solitons. The energy E and the charge Q are 
shown as functions of the frequency uj for the spherically symmetric Q-ball (a); Q-vortex (per 
unit length) (b); and Q-wall (per unit area) (c). 

The regular vortex solution exist for a frequency within the interval (j7]). It has prop- 
erties similar to the spherically symmetric Q-ball, in particular as u"^ — )■ the Q-vortex 
expands, its charge and its energy diverge. However, it is plausible that before approaching 
this limit the vortex may become unstable with respect to linear perturbations, it could 
decay into Q-balls whose energy per unit length is smaller. In the opposite limit cu^ — )■ cu^ 
the vortex also expands, however, as seen in FigJU in this limit the value of the field at 
the origin $(0) decreases and both the integrated energy per unit length 

oo 

E = 2n j dpp [{g'f + u^g^ + U{g)\ (12) 
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and the integrated charge per unit length approach constant values, as illustrated in 
Figffl (b). 

Another configuration can be constructed by imposing condition of planar invariance 
[7]. Then the field $ does not depend on p and we can use the ansatz 



$ = h{z)e 



(13) 



Substituting this ansatz into the field equation we obtain the equation without the 'friction' 
term 

which corresponds to the globally regular 'Q-wall' solution discussed in a different model 
in 13. 

(a) (b) 
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Figure 2: Energy-charge ratio versus w (a) and the energy versus charge are shown for the 
spherically symmetric Q-ball (b), Q-vortex (per unit length) (c) and Q-wall (per unit area) (d). 
A straight dashed line E = mQ indicating the margin of stability is drawn. 

For this configuration the field h{z) smoothly interpolates between some finite value on 
the z-axis h{0) and zero value as 2; — > oo. Indeed, linearisation of the equation (11 4p around 

h = gives asymptotically /i ~ e so the field decays exponentially fast. The 

solution also is regular at 2; = since h'{0) = 0. 
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As seen in FigJH the value of h{0) depends on co; as u"^ — )■ it smoothly decreases 
towards zero although the Q-wall expands. In this limit the energy and the charge of the 
Q-wall per unit area tends to zero, as illustrated in Fig[T](c). As cu^ starts to decrease the 
energy density of the wall per unit area increases forming a single maximum at the origin. 

(a) (b) (c) 



Figure 3: The energy isosurfaces of the stationary non-spinning Q-ball (a), Q-vortex (b) and 
Q-wall (c) are shown for oj = 0.8. 

In the opposite limit u"^ — the wall rapidly expands forming new vacuum ^^ac = 
1 > $(oo) around the origin, as seen from Figures O O Then both the charge Q and the 
energy E per unit area are growing without limit (Fig. [D (c)), so in all the cases there is a 
different behavior as u"^ tends to = m^. 

The fundamental solutions of the model ([T]) are stable under condition E < mQ [1]. 
In Fig. [2] we have plotted the energy E as a function of Q and the ratio E/Q versus 
uj for the different types solutions. Here we have also plotted a straight line E = mQ 
indicating the region of stability of the configurations with respect to decays into quanta 
of the scalar field, so the fundamental Q-ball becomes unstable as u approaches the critical 
value Ucr ~ 0.89. Charge and energy then approach a minimum forming a typical spike 
profile (see Fig 12], frame (b)) which indicates the saddle node bifurcations and instability 
of the upper branch. On the other hand, both the Q-vortices and Q-walls remain stable 
over most of the range of values of u (see Fig. [21 frames (c) and (d)). Note also a peculiar 
kink which appears on the upper branch of the frame (d) representing the Q-wall energy 
per unit area versus charge per unit area. 

3 Coupled Q-balls 

Clearly the relation between the energy and charge of a stationary non-topological soliton 
of Q-ball type depends on the explicit form of the potential. In the case of a fiat poten- 
tial inspired by supersymmetric models [U [13] the energy/charge relation is E ^ Q^^^, 
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Figure 4: The value at the origin of the profile functions f{r), g{p) and h{z) of the fundamental 
Q-ball, Q-vortex and Q-wall solutions, respectively are shown as a function of oj. 
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Figure 5: The profile of h{z) (a) and the energy density distribution (b) for the Q-wall solutions 
is presented for uj"^ = 1.095,0.64,0.105. 

which in 3d precisely matches the topological lower energy bound for the solitons of the 
Faddeev-Skyrme model |14|. Note that this similarity holds not only in 3 spatial dimen- 
sions because for the d-dimensional SUSY Q-balls the energy-charge relation is given by 
E ~ Q'^''*^'^'''^^ [l3]. On the other hand, for the generalized Faddeev-Skyrme model with a 
map of degree k $ : M^''"^ S"^^ C R^''^^ and the linking number Q G Z = 7r4fc_i(5'^'') 
the generalized Vakulenko-Kapitanski bound is E > Q^~^/^^'^'> = [T6]. Another in- 

teresting similarity between these models is that the spectrum of solition solutions in both 
cases includes not only the fundamental localised solitons, the Q-balls and the hopfions, 
respectively, but also extended objects like vortices [El [3] and walls [T8l 17]. 

Recently Radu and Volkov [Ij noted that there is an interesting possibility to establish 
a link between some generalization of the model ([1]) and soliton solutions of the Faddeev- 
Skyrme model in 3 spacial dimensions. They suggested to consider the generalized spinning 
ansatz for the field $(r) which also includes an independent phase \l/(r, ^): 

$ = /(r, oy-'t-m^ir,e)+nv> ^ ^x{r, 6) + iY{r, ^)]e^'^*+"^, n, m G Z. (15) 
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Figure 6: The energy density of the Q-wall configuration is shown for uP' = 0.64 (a) and = 
0.105 as function of x^z coordinates. 



Similar to the parametrization of the first two components of the axiaUy symmetric hopfion 
field [17J, the phase function \E'(r, ff) increases by l-n after one revolution around the closed 
contour which consists of the 2;-axis and a semi-circle whose radius expands to infinity. The 
phase + then increases by 27rm times as one moves along this contour and by 2?™ 
times as one moves on a circle around the ^-axis. Such a 'twisted' Q-ball is characterised 
by 2 integer winding numbers n and m, which, however do not have a topological meaning. 

Thus with our choice of the parameters of the potential, the fields X(r, &) and F (r, &) 
of the twisted Q-ball are coupled via the interaction lagrangian 

L,„, = 3X2r2(x2 + r2_4/3) (IQ) 
and the energy density of the stationary system in cylindrical coordinates read 



+ Ui[X] + U2[Y] + L 



int 1 



where 

Ui[X]=aX^ + bX'^ + cX^; U2[Y] = aY^ + bY^ + cY^ . (17) 

Actually we are considering two coupled configurations of the Q-ball type with two similar 
sextic potentials ([2]). This coupling, however, is rather restrictive, unlike in the previous 
examples considered in [20] and in the parameters of the potentials f|T7|) and the 
coupling constant are fixed by the structure of the model. Consequently the parameter u) 
and the winding n are the same in both sectors, so the U{1) charge of the configuration is 

Q = 2co [ d^x{X^ + Y^) . (18) 
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3.1 Coupled non-spinning n = solitons 

The ansatz f|T5l) provides a natural way to couple two fundamental non-topological solitons. 
First, let us consider coupling of two fundamental solitons of the Q-ball type provided by the 
ansatz f|T5l) with winding number n = 0. The Euler-Lagrange equations arising then from 
the variations of ([1]) with respect to the functions X{p, z) and Y[p, z) have been integrated 
by numerically imposing the boundary conditions, which respect finite mass-energy and 
finite energy density conditions as well as regularity and symmetry requirements. 

The numerical calculations are performed employing the package FIDISOL/CADSOL, 
based on the Newton- Raphson iterative procedure [19] . We solve the system of two coupled 
nonlinear partial differential equations numerically, on a non-equidistant grid in p and z, 
employing the compact coordinates x = p/(l -|- p) G [0 : 1] and ?/ = 2;/(l -|- 2) G [0 : 1]. 
Typical grids used have sizes 85 x 70. The relative errors of the solutions are of order of 
10~^ or smaller. 

It turns out that the corresponding non-spinning = coupled two-component configu- 
rations do not exist for the constituents of different geometry, e.g., we do not find solutions 
which would represent the Q-wall coupled to the Q-vortex/Q-ball. We do however find 
two-component coupled solutions of the same type with z) = Y{p, z) interacting via 
the potential term (|T6|) . We could expect the solutions may exist only in the parameter 
range < < ([7]) where = 0.1 and uj\ = m? = 1.1. 

Considering two spherically symmetric coupled Q-balls we found that the corresponding 
integrated energy and the U{1) charge Q of the system depend on u in almost exactly the 
same way as in the case of the single component model considered in the previous section. 
Furthermore, the energy is equally distributed between two components, so the 2-Q-ball 
system remains spherically symmetric although the behavior of the profile function /(O) at 
the origin is different. This is illustrated in Fig. [TJ frame (a), where we exhibit the energy- 
charge ratios E/Q of the 2-Q-ball system and the single spherically symmetric Q-ball as 
functions of u. Fig. [TJ frame (b) displays the evolution of the profile function /(O) at the 
origin for these systems. 

In the limiting cases u — t- w+ and u ^ U- behavior of the spherically symmetric 2-Q- 
ball system is similar to the case of the single spherically symmetric Q-ball. However the 
potential of interaction of the components ( !T6|) has a different structure. For 0.90 < u < 
it is attractive and it has a minimum at the origin, as seen in Fig. [H] (b). The depth of the 
spherical potential well increases as u decreases. For u < 0.90 the potential of interaction 
is repulsive at the region around r = 0, it has a minimum at r = ro 7^ 0. Thus the energy 
density distribution of the 2-Q-ball coupled system is forming a shell structure (see Fig. |8] 
(a) where we presented the energy isosurface at a; = 0.50). 

For two axially symmetric non-spinning coupled Q-vortices and Q-walls we observe a 
similar pattern. An initial configuration which satisfies the proper boundary conditions, 
rapidly converges to the two-component solution with X{p) = Y{p), however the energy 
and the charge per unit length (area) of the coupled system almost coincide with corre- 
sponding values of the energy and the charge of the single- component fundamental solution 
over all range of values of u presented in Figs. HJ [21 Evidently, the energy/charge relation 
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Figure 7: Energy-charge ratio versus u is shown for the coupled spherically symmetric 2-Q-ball 
system and for the single Q-ball (a). A straight dashed line E = rriQ indicating the margin of 
stability is drawn. The value of the profile functions /(O) at the origin is shown as function of co 
(b). 

on the upper branch is rather different from that of the solitons in the Faddeev-Skyrme 
model. Thus, the energy and the charge are equally distributed between the components 
and the values of the profile functions on the symmetry axis (symmetry plane) is smaller 
than they are in the case of the single-component Q- vortex/ Q- wall (see Fig. [9] (a,b)) 

Again, as w — t- a;+ and a; — t- the behavior of the coupled axially symmetric 2- 
Q-vortex system and 2-Q-wall configuration follow the familiar pattern of the single Q- 
vortex/Q-wall respectively. As u decreases, both components of the coupled system ex- 
pand, the potential of interaction then is repulsive on the symmetry axis/symmetry plane 
and it has a minimum at p = po(i^) and z = zq, respectively. Thus the energy density dis- 
tribution has the shell structure (see Fig. [10] (a)). Then both components rapidly expand 
forming a new vacuum around the 2;-axis/x — y plane, as seen from Fig. [10] (b) (cf similar 
behavior of the Q-wall, Figs. \5\ Due to severe numerical difficulties encountered here, 
we could not clarify the properties of the critical solution as u approaches a;_, in particular 
it is unclear whether the bounded 2-Q-vortex system remains stable in this region or it 
decays into Q-balls with lower energy per unit length. 

3.2 Coupled spinning n = 1 solitons 

The ansatz f[T5|) also describes spinning configurations with n ^ 0. Indeed, it was demon- 
strated in [3| that the single component model admits spinning, axially symmetric 
generalizations for the spherically symmetric Q-balls. The properties of these solutions 
were analysed in [201 [SI [SI [21] ■ The energy density of these spinning Q-balls is of toroidal 
shape and they possess the angular momentum J = nQ. It has also been shown by Volkov 
and Wohnert [3j that for a given n both parity-even (i.e. symmetric, P = +1) and par- 
ity odd (i.e., antisymmetric, P = —1) solutions exist. Furthermore, there are radial and 
angular unstable excitations of the fundamental Q-balls which are related to the spherical 
harmonics [31 12D]- The radial excitations are parametrised by the number of nodes of the 
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(b) 
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Figure 8: The energy isosurface of the 2-Q-ball system is shown for uj 
of interaction is plotted as function of r for w = 1.04 and uj = 0.50 (b). 



0.55 (a). The potential 



scalar field G Z. Also the spinning excited Q- vortex solutions are known [3]. 

However, we observe that no spinning Q-wall configuration with n 7^ is likely to exist. 
This agrees with the physical intuition based on a symmetry argument. However we have 
found some evidence that there are families of generalized Q-wall solutions spinning around 
the p axis. This study will be reported elsewhere. 

A general property of the spinning 2-soliton solutions parametrised by the ansatz f|T5|l 
for the system ([1]) is that the interaction between the components rapidly drives the system 
to the degenerated state X{p,z) = Y{p,z). We restrict our discussion to two particular 
cases considering first the system of parity-even X-component in the sector with n = 
l,k = and secondly, the parity-even angularly excited F-component with n = 1, k = 0. 
This system would provide an example of a " non- twisted" (1,0) Q-ball (|T5i) with winding 
number m = 0. Similarly, one might try to find a "twisted" Q-ball solution, employing the 
parity-even X-component in the sector with n = l,k = 1 and the parity-odd angularly- 
excited F-component [1] . These configurations serve as a first guess to obtain solutions of 
the Euler-Lagrange equations. 

The numerical calculations clearly show that it is energetically favorable in both of these 
cases to have two component Q-balls sitting on top of each other with both constituents 
having identical geometry. Thus, the system converges to the two copies of the rescaled 
axially symmetric angularly-radially excited Q-ball which are parity-even and asymmetric, 
respectively (see Fig. [TT]) . 

As seen in Fig. [TT](a) for the parity even axially symmetric spinning (1,0) 2-Q-ball con- 
figuration the fields have three almost identical maxima at a finite distance from the ^-axis 
with one maximum in the equatorial plane and two other maxima located symmetrically 
with respect to the x — y plane. The energy-momentum and charge densities then exhibit 
a system of three tori presented in Fig. [T21 left frame. 
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Figure 9: 2-Q-vortices (a) and 2-Q-wall (b) coupled system: The value of the profile functions 
g{0) and h{0) at the origin are shown as function of uj together with the corresponding function 
of the single component model. 




Figure 10: The energy isosurface of the 2-Q-vortex system at Tg = 0.21 (a) and the energy 
density of the 2-Q-vortex configuration as function of x, z coordinates (b) are shown for uj = 0.43. 

For the asymmetric spinning (1,1) 2-Q-ball configuration the functions X(p, 2;) = 
z) have a bit more complicated stricture which can be interpreted as superimposed 
field of two excited Q-balls. For the first radially excited spinning symmetric Q-ball, the 
amplitude has one node in the equatorial plane and the energy and charge densities form 
a system of two concentric tori, for the second angularly excited antisymmetric spinning 
Q-ball, the amplitude also has a node in the equatorial plane and the energy and charge 
densities show two maxima located symmetrically with respect to the x — y plane as shown 
in Fig. [T2I right frame. 

In Fig. [12] we plotted the energy E of the (1,0) and (1,1) composite 2-Q-ball configu- 
rations as a function of Q (frame (a)) and the ratio EjQ versus lo (frame (b)). As before, 
the straight line E = niQ indicates the region of stability of the configurations with re- 
spect to decays into quanta of the scalar field. In both cases the behavior of the 2-Q-ball 
systems is almost identical to what we observed for the single axially symmetric angularly- 
radially excited Q-balls. Thus the (1,0) 2-Q-ball becomes unstable as u approaches the 
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Figure 11: The profile of tiie function X[p,z) = Y{p,z) is sliown for the parity even "non- 
twisted" (1,0) Q-ball (a) and for the mixed parity "twisted" (1,1) Q-ball (b) for a; = 0.70. 

critical value Wcr'* ~ 0.835 and the (1,1) 2-Q-ball is unstable with respect to decays as 
u > ooil^ ~ 0.847. Both values are smaller than the critical frequency of the fundamental 
spherically symmetric Q-ball. 

Charge and energy then approach a minimum forming a typical spike profile (see FigJT3t 
frame (b)) which indicates the saddle node bifurcations and instability of the upper branch. 
On the other hand, the composite spinning 2-Q-balls parametrized by the ansatz fllSp 
are not stable with respect to decay into the fundamental spinning Q-balls. Numerical 
calculations reveal this instability as u decreases below a certain threshold value which is 
much higher than the cu^ = 0.1, for example the (1,0) and (1,1) 2-Q-balls are unstable 
with respect to decay into the coupled system of 2 radially excited n = 1 Q-balls as u 

(2) 

approaches the second critical value oocr ~ 0.60. 

Let us finally briefly discuss one more possibility to construct a composite system of 
two coupled non-topological solitons [201 S] • We follow closely the approach described in 
|20] . The Lagrangian density for this model is given by 

L = L(<l)i) + L($2) - A|<l>in$2p (19) 

where two copies of the Lagrangian ([1]) L($i) = — [/(I'l'il), i = 1,2 are coupled 

through the minimal interaction term. This coupling is much less restrictive than the 
example ffTSjl above. Indeed, both the coupling constant A, the frequencies uji,uj2 and the 
windings ni,?7,2 are free parameters of the model. As compared to the parametrization 
( IT^ . we expect the existence of a much richer set of possible composite solutions for the 
model ([in]). Some of them were considered in [201 S], here we would like to note that it 
is now possible, within this model, to couple two solitons having different geometry. In 
Fig. [in we exhibit the energy energy isosurfaces of some two-component configurations at 
A = 0.1 and ui = U2 = 0.90: 

(a) Q-wall { rii = 0) and the fundamental Q-ball (n2 = 0); 

(b) Q-vortex {rii = 0) and Q-wall (^2 = 0); 
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Figure 12: Isosurfaces of constant energy density of the " non- twisted" (1,0) Q-ball with Tq = 0.32 
(left frame) and the "twisted" (1,1) Q-ball with = 0.35 (right frame) are shown for oj = 0.90. 

(c) Q-vortex {rii = 0) and spinning Q-ball (^2 = 1) ('hoop' solution [1]); 

(d) Q-wall (^1 = 0) and spinning Q-ball (^2 = 2); 

(e) Q-vortex {rii = 0) and the even parity angularly excited spinning Q-ball (n2 = 
2) (double hoop configuration); 

(f) Q-wall {rii = 0) and the odd parity angularly excited spinning Q-ball (^2 = 1). 

Obviously the stability of these configurations depends on the coupling A. The inves- 
tigation of these and other two-component solutions of the model ( 1T9]) and their eventual 
bifurcations is currently underway. 

4 Conclusions 

Motivated by the recent interest in non-topological stationary solitons possessing vari- 
ous types of symmetry, we investigated properties of the corresponding solutions to the 
same simple model with sextic potential: configurations with spherical symmetry (Q-balls), 
cylindrical symmetry (Q- vortices) and planar symmetry (Q-walls). 

The numerical studies of these three types of fundamental solutions show different 
behavior of charge and energy with changing u for each type. Their dependence on u 
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Figure 13: Energy-charge ratio versus w (a) and the energy versus charge (b) are shown for the 
(1,0) and (1,1) 2-Q-balls. A straight dashed line E = mQ indicating the margin of stability is 
drawn. 

is analysed numerically in some detail. By and large this is qualitatively very similar to 
that for the spherically symmetric solutions [IJ. In all cases there is a different behavior 
as oj tends to 0;+ or a; — )• We have also found numerical evidence of an instability of 
the Q- vortices and the Q-walls in the limit where a; — )• a;_. Our numerical investigations 
indicate fairly clearly that there is no spinning generalization of the Q-wall solution, which 
would be constructed by a direct analogy with spinning Q-balls and Q-vortices. It does 
not exclude another possibility that there are generalized Q-wall solutions spinning around 
an axis in the x-y plane. 

Also, we have studied a coupled two-component system using the generalized field 
ansatz which includes an independent phase [4j and have constructed several examples 
which represent both spinning and non-spinning configurations. Previous research has 
suggested that such a system may support the existence of 'twisted' Q-balls which would 
have a certain similarity with twisted loops in the Faddeev-Skyrme model. As it turns 
out these coupled configurations do not exist if the constituents possess different geometry. 
If the constituents are of the same geometry, however, then it is energetically favorable 
to have two component sitting on top of each other. Thus, the 'twisted' system rapidly 
converges to the two rescaled copies of the single component model with the energy and 
the charge equally distributed between the components. 

However such a composite system is unstable not only with respect to radiation of the 
scalar quanta as a; — )■ a;+, but also with respect to decay into the fundamental solitons. 
Numerical calculations reveal this instability as u decreases below certain threshold value 
which is much higher than the U-. 

Another possibility, discussed in [20], is not so restrictive. A minimal interaction be- 
tween two components allows us to construct a plethora of coupled Q-balls with various 
geometry. The numerical work involved in the construction of these 2-Q-ball solutions is, 
however, a considerably challenging task. 
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(a) (b) 




Figure 14: The energy isosurfaces of the coupled system with component fields of different 
geometry are shown for ui = 0J2 = 0.90 and coupling A = 0.1. 
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